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We study the sensitivity of Tc and the pseudogap onset temperature , T*, to low fields, H, using a BCS-based 
approach extended to arbitrary coupling. We find that T* and Tc, which are of the same superconducting origin, 
have very different H dependences. This is due to the pseudogap, Apg, which is present at the latter, but not 
former temperature. Our results for the coherence length ^ fit well with existing experiments. We predict that 
very near the insulator ^ will rapidly increase. 

PACS numbers; 74.25.Ha, 74.60.Ec, 74.25.-q COnd-mat/01 0361 4 



One of the central questions in understanding the under- 
doped cuprates is the extent to which the superconducting 
phase is described by BCS theory. Recent experimentsBij in- 
dicate that the pseudogap persists below Tc in the underlying 
normal density of states. Thus, the fermionic excitation gap 
A is to be distinguished from the order parameter A^c- These 
two energy gaps mirror a distinction between the two temper- 
atures T*(the pseudogap onset), and Tc(the superconducting 
transition), which behave differently as a functinn of hole con- 
centration X as well as of magnetic field, _HaUu. Indeed, T* 
and Tc are, respectively, weakly and strongly dependent on 
H in the well-estabUshed pseudogap regime. Moreover, the 
distinction between tbfS^ temperature and energy scales has 
been frequently citedQ aEl as evidence that they have different 
physical origins. 

In this paper we provide a counter argument to this widely 
stated inference by demonstrating that these crucial magnetic 
field effects in the pseudogap phase, are entirely compatible 
with superconductivity as origin for both T* and Tc- Our 
approach is based on an extended version of BCS theory, in 
which the attractive coupling g is contemplated to be strong 
enough so that pairs begin to form at a higher-temperature 
T* than Jhe Tc at which they Bose condenseoB. We have 
shownETtil that as a necessary jCpmpefluence A ^ A^c. More- 
over, our work has emphasizedBE3£j that a (pseudo)gap in the 
fermionic spectrum at Tc is deleterious for superconductivity. 
Thus, as observed experimentally, as a function of decreasing 
X, Tc decreases as the pseudogap or T* grows. 

A calculation of the field dependence of Tc {i.e., Hc2) is an 
important problem in its own right, (i) This is the only way to 
provide a precise interpretation of the "coherence length" i^, 
which Jiie demonstrate here is very different from that of BCS 
theorytJ. (ii) An analysis of Hc2 is tantamount to arriving 
at a reformulation of the microscopically deduced Ginzburg- 
Landau (GL) free energy up to quadratic terms, which must 
necessarily incorporate the presence of a non-zero pseudo- 
gap at (iii) Because the "competing order parameter" 
scenarictl'E£l also addresses the observation that A ^ Asc, as 
well as the competing x dependences of Tc and T*, magnetic 
field effects may provide a unique testing ground for distin- 
guishing between these two scenarios. 



Indeed there is a rather close similarity in the structure of 
our zero-field theory to the phenomenplogy deduced from 
thermodynamical data by Loram et aO. Our mean field 
calculationsLiU show that the gap equation for T < Tc re- 
duces to the usual BCS form, but with a new quasi-particle 
dispersion 
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Here ipk is associated with the pairing symmetry. It follows 
from this that the larger is Apg{Tc), the lower the transition 
temperature Tc- In contrast to the work in Refs. |6| and ^s], 
here Apg{T), is determined self-consistently and derives from 
the presence of a strong pairing attraction. Moreover, this 
pseudogap, Apg, persists below Tc. 

The conclusions of this paper are relatively straightforward 
and we begin with a simple intuitive argument to address iJc2- 
Consider the Ginzburg-Landau free energy functional near Tc 
to quadratic order in Age, in a finite field 
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Here tq describes how the system approaches the critical point 
with varying temperature, and rj^ is the stiffness against spa- 
tial variations of the order parameter The mean-field behavior 
of To near Tc yields to{T) ^ fo(l - T/Tc). It follows that 
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where ^ is the zero temperature coherence length. A rough 
extrapolation yields -ffc2(0) ~ $o/(27r^^). In the small g 
(i.e., BCS) case tq = N{0), the density of states per spin at 
the Fermi surface, and r/^ = N {0)7 ( (3) / ASn^ (vp /Tc)^ ■ The 
squared coherence length ^^^5 — 7C(3)/487r^(t;_F/Tc)^ is 
determined by the stiffness 77^ with tq cancelUng the density 
of states. 

In contrast, in the strong coupling case the pseudogap mod- 
ifies the fermionic quasiparticle dispersion through a replace- 
ment of ek by i?k and thereby suppresses fo. Moreover, the 
stiffness 77^ (which is relatively insensitive to the energy scale 
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of the pseudogap), decreases due to the diminishing pair size. 
There is, thus, a competition between the numerator and de- 
nominator in Eq. (^. The decrease in fo dominates at suffi- 
ciently large coupling g, resulting in an extended flat region 
followed by an eventual growth of £^ with increasing g. The 
latter reflects the-approach to the ideal "boson" limit, where 
Tc is suppressecL-3 to zero at any H ^ 0. 

Next, we provide a microscopic derivation of the param- 
eters in Eq. (^j) (and a related counterpart for T*). A cen- 
tral theme in our paper is that both the field sensitivity of 
Tc and of T* in small H can be studied through the zero- 
field normal state pair propagator, or the inverse t-matrix 
t~^{Q) — l/ff + x(Q)' where Q — (q, ^i) is a four-vector. 
The coefficients in Eq. (||) will be shown to arise from an ex- 
pansion of t^^ in the momentum components perpendicular 
to the field (indexed hy i, j — x,y) as 
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To = ^ + X(0, 0), = 2 ydet [d,,d,^x{Q)] 
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where we have generalized from Eq. (^ to include possible 
anisotropy. 

tion for our choice of x{Q) is essential, 
this paper combines a Green's function 
with a generalization of the BCS ground 
This formalism allows for Fermi- and 
non-Fermi-liquid (i.e., Apg{Tc) / 0) based superconductiv- 
ity, according to the size of g, with self consistently deter- 
mined chemical potential /i. In the present paper all technical 
issues of this decoupling scheme can be simply by-passed, and 
the results obtained are not only intuitive, but rather general. 
All that is needed here is the observation that the pair suscep- 
tibility x(q,f^) = X(Q) = J:kG{K)Go{Q - K)^i_ 
Thus 
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x(q,o) = X!'^k-q/2 = 
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Here and are the usual BCS coherence factors, and 



= (1 + {k/koYr 



or 



{cos{kxa) — cos{kya)) , for 



wave pairing in 3D jellium or d-wave pairing on a quasi- 2D 
lattice, respectively. That there is one full Green's function 
(G) along with one bare Green's function (Go), reflects the 
structure of the BCS gap equation, which introducesEJ this 
xiQ) form, (with integrand proportional to the usual Gor'kov 
F function). All numerical calculations in this paper are based 
on Eqs. (|), and (|), given Apg = Apg{Tc). Although 
here we proceed more self consistently, our analytical scheme 
for computing the various energy scalesEJ can be by-passed, 
if Apg and Tc are pre-determined, e.g., fitted to cuprate exper- 
iments. 

We next turn to T*, where the Fermi liquid begins tp-hr^ak 
down; this is associated with the onset of a resonanceE3E 
t{Q), as g becomes sufficiently large. Detailed numericsE 
based on the coupled Green's function equations show that to 
a good first order approximation this resonance temperature 



can be deduced from the condition l/g + xo(0, 0) = 0, where 
Xo is given by Eq. (^ with substituted for i?k- Indeed, quite 
generally, at T > T*, the ^-matrix can be well approximated 
by using xo{Q) in place of x{Q)- 

Magnetic field effects can be readily included into our for- 
malism. We begin with a derivation of T*{H) to linear order 
in H. Our Hamiltonian consists of the field dependent kinetic 
energy term along with the usual two-body pairing interaction 
of strength Vk_k' — SV-kVli' where we now include interac- 
tions between pairs of non-zero net momentum. Fluctuations 
of pairs with finite momentum (q, k) are characterized by the 
correlation function 

D{q,k;q,k') = {T b{q, k,T)b'^ {q ,k' ,t')) , 

where fc, r) = e(«-^^)-c,/2_feC,/2+fee-(«-'^^)^ 

Summing ladder diagrams leads to a Dyson equation D = 
GqGq — GqGqVD with solution 
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where xo = G^Go is the counterpart pair susceptibility 
for iJ / 0, and the field-dependence of the bare electron 
propagator GqOS implemented using the semiclassical phase 
approximationElelevating both xo ™d J^k k' E'^k'-Pk' to in- 
tegral operatorsE^ whose eigenvalues satisfy Eq. (^. This 
approximation allows the calculation of the eigenvalues of 
Xo from the zero-field pair susceptibility xo in the regime 
T ^ eH/ mc. The field-induced relative phase shift between 
electrons in a pair renormalizes the interaction Vk.ks but the 
effect is quadratic in H and is therefore ignored here. The 
pairing resonance temperature T* is defined by the appear- 
ance of an eigenvalue Ho = —.9^^ of xo^ which causes Eq. (||) 
to diverge. We define parameters rf^ and Tq — f^{l~ T/T*) 
analogous to those which appear in Eq. (0) and obtain 
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IIo = Tn 



*2 . ^i? = 

c 
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which defines T*{H) to linear order in H. 
T*{H) atff = Ois 



1 dT* 
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which leads to the associated "coherence length" = 
■q*"^ /tq. The stiffness if ^ can be explicitly evaluated using 
the zero-field pair susceptibility, and 
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T=T' 



Here we have included a contribution from the temperature 
dependence of /i. In the weak coupling case, the chemical 
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FIG. 1. Coherence lengths ^ (for Tc, solid) and ^* (for 
T* , dashed) with variable coupling g for s-wave jellium; here 
Qc = — 47r/mfco and units are fc^^. See Eqs. (jsj) and (j^. Inset 
shows the zero field behavior of T* , Tc, and Apg from RefO, where 
T and A are in units of Ep. 




FIG. 2. Coherence lengths at Tc (solid) and T* (dashed) versus 
variable coupling g for a d-wave lattice at density n = 0.85. ^ is4fi 
units of a. The inset plots the zero field energy scales from Ref Ji3, 
and T, A are in units of 4t|i . 



potential is pinned at Ep, and we recover the (s-wave) BCS 
limit To* = - /'(ek)¥'L« « 7V(0). 

The dashed line in Fig. ms a plot of the slope of T* with 
H, as a function of the coupling g, for the case of s-wave 
jellium. This should be compared with the inset, from Ref. |l3[ 
which illustrates the suppression of Tc by Apg for all g in the 
fermionic regime; (ultimately, when /i becomes negative, Tc 
starts to increase again). When pseudogap effects are weak, 
T* is essentially the same as Tc and depends strongly on H. A 
stronger pairing interaction g causes Apg to increase and the 
pair size to decrease; the latter effect diminishes the stiffness 
and causes T* to be weakly field dependentEJ. 

We turn next to Tc{HX-&Vid note that a solution of the cou- 
pled equations of motionllS (as was done in the zero field case) 
appears prohibitively difficult. Nevertheless, based on the 
above observations that (i) in_zero field T* scales cather well 
with Apg (both theoreticallyllj and experimentallyH,) and (ii) 
that T* is very weakly field dependent in the well-established 
pseudogap regime, we infer that Apg = Apg(Tc) is weakly 
H dependent. This assumption, along with the semi-classical 
phase approximation for the full Green's function, are the only 
essential assumptions made here. The weak ii dependence in 
Apg appears compatible with experimentQ' 



Od and underlies 
a GL formulation (Eq. (Q)) in which only the superconduct- 
ing order parameter is coupled to the magnetic field. In this 
way, for the purposes of computing Hc2, the pseudogap enters 
as a relatively rigid band structure effect, which is accounted 
for by introducing the full pair susceptibility into the standard 
Hc2 formalismcj. This approach necessarily yields the correct 
H ^0 result for Tc. 

It also leads naturally to Eq. (Q), from which one deduces a 
rather complicated expression (not shown) for rj^ along with 



(10) 



where we have omitted contributions from the temperature de 
pendence of /i and Apg near Tc- Eq 



(10) contains the es- 
sential physics introduced by the pseudogap. The summa- 
tion essentially measures the E ^ density of states which 



is depleted by the pseudogap at strong coupling, leading to 
To g^-^va/Ta^ Moreover, we have shown analytically that 
the neglected terms further suppress fo, and, therefore, do not 
qualitatively change our results. 

Figures [l| and ^ indicate how the characteristic low-field 
slopes behave as a function of coupling g, for the s-wave (jel- 
lium) and d-wave (lattice) cases, respectively. For the latter, 
the bare band dispersion was presumed to be ek — 2i|| (2 — 
cos(fca;a) — cos(fcj,a)) — 2tj^(l — cos(fc_Ld)) — where a is 
the lattice constant in the plane and d is the distance between 
layers. The insets in each figure summarize the|Zero field re- 
sults for Tc, r*, and Apg, calculated elsewherell3. It should 
be noted from the insets that there is an extended regime of 
coupling constants over which BCS behavior (Apg = 0, and 
Tc — T*) is obtained. However, beyond a "critical" coupling 
(corresponding roughly to where bound states of the isolated 
pair occur,) the pseudogap becomes non-zero, and Tc is dif- 
ferentiated from T*. 

It is clear from the figures that the slopes of T* and Tc 
are identical at weak coupling and become progressively more 
distinct as the coupling is increased. The two associated stiff- 
ness parameters decrease with coupling in a similar way, but 
the field dependence of Tc is enhanced by the strong suppres- 
sion of To relative to Tq . The competition between the nu- 
merator and denominator in Eq. (^ leads to a length scale 
which is relatively constant over a rather wide range of moder- 
ate g; then at sufficiently strong coupling begins to increase 
due to the reduction in the density of states (fo) associated 
with the growth of the pseudogap. 

The key role of fo in determining the squared-coherence 
length highlights the fact that at strong coupling, the vari- 
ous length scales of the system must be carefully distinguished 
even though they are identical in BCS theory. This is illus- 
trated in Figs. B(a) and (b) for the s- and d-wave cases, re- 
spectively, in which we compare normalized stiffness, fj = 
with the BCS coherence length at Tc, and ^, the 
calculated coherence length at Tc- The density of states ef- 
fect is evident in both the s- and d-wave cases, as the upturn 
of ^ at strong coupling contrasts sharply with the tapering of 
the stiffness fj. Note that ^bcs exhibits unusual structure in 
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FIG. 3. Coupling constant dependence of the key length scales: 
the calculated coherence length ^ obtained from Hc2 (solid), f) 
(the normalized stiffness; dashed), and the BCS coherence length 
(oc vf/Tc\ dot-dashed), for the s-wave jellium case (a) and the 
d-wave lattice case (b). 
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hole concentration (x) 
FIG. 4. Magnetic length scales associated with Tc and T* as a 
funotipn of hole concentration in the cuprates. Zero field results from 
RefJij are shown in the inset ( Apg is a factor of 2 smaller than 
experiment due to our convention for ifk )■ 



Fig. |(a) which arises from the non-monotonic behavior of Tc 
evident in the inset of Fig. [l|. In Fig. ||(b), the divergence of 
£,BCS at large g is "accidental" and derives,from the vanishing 
of Tc which arises from pair-localizationEj (see Fig. ^, inset). 
BCS relations, therefore, are very misleading when one tries 
to infer the behavior of ^ at strong coupling from other length 
scales in the system. 

In order to map the coupling g onto hole concentration x, 
we introduce an x-dependent hopping roatrix element t\\ [x] = 
Iqx associated with the Mott transitiontj. We presume in the 
absence of more microscopic information that 5 is x indepen- 
dent, leaving one free parameter in our theory g/to, chosen to 
optimize a fit to the phase diagram (Fig. ^ inset.) One could, 
alternatively, bypass all assumptions concerning Mott insula- 
tor physics, if one instead used experimental data in place of 
the calculations in the inset. The resulting behavior of ^ and ^* 
(Fig. ^) would be essentially the same. Over most of the range 
of X, ^ is relatively constant, as seems to be observed experi- 
mentally, ancLits magnitude is within a factor of two or three 
of experimentEB. As in expetdment, T* is found to be less field 
sensitive in the underdopecm than overdopeJB regimes. As the 
insulator is approached, ^ rapidly increasesE3 while ^* contin- 
ues to decrease. We have thus demonstrated that the different 
observed field dependences of T* and Tc (for both under- and 
over-doped cuprates) are contained in our theory, in which the 
pseudogap is associated with precursor superconductivity. 
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Note Added. — After this work was completed weJearned 
of related experimental studies by Shibauchi et al.Lj which 
show strong similarities to our theoretical predictions. 
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